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ABSTRACT 



in 

> 

■ Context. A significant number of the known multiple exoplanetary systems are containing a pair of giant planets engaged in a low order mean 
i^-j motion resonance. Such a resonant condition protects the dynamics of these planets resulting in very stable orbits. According to recent studies 
£f~^ . the capture into a resonance is the result of a planetary migration process induced by the interaction of the planets with a protoplanetary disk. 

If the migration is slow enough (adiabatic) next to a mean motion resonance, the two planets will also be in apsidal corotation. 
\Q) Aims. The recently refined orbital parameters of the system HD 128311 suggest that the two giant planets are in a 2:1 mean motion resonance, 
f*"^ ■ however without exhibiting apsidal corotation. Thus the evolution of this system can not be described by an adiabatic migration process alone. 

We present possible evolution scenarios of this system combining migration processes and sudden perturbations. 
^ , Methods. We model migration scenarios through numerical integration of the gravitational N-body problem with additional non-conservative 
I ■ forces. Planet-planet scattering has been investigated by N-body simulations. 
O Results. We show that the present dynamical state of the system HD 1283 1 1 may be explained by such evolutionary processes. 
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1. Introduction J2002l) : lBeauge et alJ J2006). in the case of a sufficiently slow 

migration process the resonant planets should also be in apsidal 
Among the 19 multi-planet systems found to date about a third corotation . In order t0 resolve the above d i scre pancy we con- 

struct mixed evolutionary scenarios of migrating planetary sys- 



are e ngaged in a low order mean motion resonance JVoet et alJ 



|200J>. The most prominent case is the exact 2: 1 resonance of tems incorporating m i gration and othe r additional perturbative 

the two outer planets in GJ 876. There the orbital elements are effects For Qur invest igation we use N-body num erical inte- 

very well determined, due to the short periods of the planets of ^i ons containing a i so non-conservative forces llLee &Pealel 

* 30 and 60 days ( |Laughlin et alj2005|) . The formation of reso- gjjj^ which haye be e n tested on full h ydrodynamic evolutions 

nant configurations between planets must be due to dissipative of embedded planets faev et all2004 . In this letter we report 

processes altering the semi-major axis. For the system GJ 876 our findings in modeling the behavior of the reS onant exoplan- 

the orbits are in apsidal corotation and both resonant angles etarv system HD 128311 
librate with small amplitudes, a condition which can best ex- 
plained by a sufficiently slow and long lasting differential mi- 

gration process induced by the interaction of the planets with a 2 Qrbjta| so | ution and jts stab i| ity 

protoplanetary disk ( Sne llgrove et al.l200U Lee & Pealel2002l 

iKlev et al1l2005l) . Hence, the occurrence of resonances consti- It is mentioned by IVo gt et all J2005h that the best fit to the 

tutes a very strong indication of migration in young planetary radial velocity curve of HD 128311, based only on the use 

systems, in addition to the hot Jupiter cases. of Keplerian ellipses, leads to orbital data resulting in unsta- 

Recent analysis of the planetary system HD 128311 sug- ble behavior. Therefore they present also an alternate fit for 

gests that two giant planets are engulfed in a 2: 1 mean mo- HD 1283 1 1 calculated by using three-body gravitational inter- 

tion resonance, h owever without exhibiting apsid al corotation actions. The orbital data determined by this alternate fit (listed 

JVogt et all2005l) . On the other hand, according to lLee & Pealel in Table 1) result in stable orbits, where the two giants are in 

a protecting 2:1 mean motion resonance (MMR). We perform 

Send offprint requests to: Zs. Sandor three-body numerical integrations using these initial conditions 
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Table 1. Orbital data of HD 12831 1 as provided bv lVogt et al] 
(2005). M denotes the mean anomaly. 



planet 


mass [My] 


a [AU] 


e 


M [deg] 


m [deg] 


b 


1.56 


1.109 


0.38 


257.6 


80.1 


c 


3.08 


1.735 


0.21 


166.0 


21.6 



of lVoPt etalJ ll2005). as listed in Table 1 . Only the resonant an- 
gle ©i = 2A 2 - Ai — vj\ librates around 0° with an amplitude 
of ~ 60° (Fig.[T] top), while both 2 = 2A 2 - Ai — vj 2 and 
Am — w 2 — tu\ circulate. Here Aj are mean longitudes and nr, 
are longitudes of periapse, both numbered from the inside out. 
This implies that the systems, although engaged in a MMR, 
is not in apsidal corotation. Moreover, the eccentricities show 
also large oscillations (Fig.[T] bot tom). 

As described by IVogt et alJ ( 120051) . due to the relatively 
high stellar jitter, there exist a large variety of dynamically dis- 
tinct stable orbital solutions resulting in the current radial ve- 
locity curve. The authors present the results obtained by a self- 
consistent two-planet model (taking into account the mutual 
gravitational interactions) and a Monte-Carlo procedure. As fi- 
nal outcome, they have found that only very few samples ex- 
hibit stable behavior without resonance. The remaining stable 
solutions are in 2: 1 MMR, in the majority of the cases (~ 60%) 
only ©i librates while ©2 circulates, in the remaining cases ©2 
librates as well. Hence, we assume that their derived parame- 
ters in Table 1 are an accurate representation of the dynamics 
of the system. 

3. Evolution through a migration process 

3.1. Adiabatic migration 

As shown in the previous section, the planets in the sys- 
tem HD 128311 are presently engaged in a stable 2:1 MMR. 
Although the actual orbital parameters do not exhibit apsidal 
corotation, the existence of the resonance suggests that the sys- 
tem has in the past gone through a migration process. 

The migration of a single planet can be characterized by the 
migration rate a/a and the eccentricity damping rate e/e. Here, 
we use the corresponding e-folding times for the semi-major 
axes and eccentricities: r fl and r e , respectively. The relation be- 
tween the damping rates and e-folding times is a/a = -l/r fl , 
and simi larly for the eccen tricities. Investigating the system 
GJ 876. iLee & Peald J2002I) have found that for a sufficiently 
slow migration the final state of the system depends only on the 
ratio K of the e-folding times K - T a /r e . 

From hydrodynamical calculations we know that the order 
of magnitude of K is close to unity and it reflects the phys- 
ical properties (i.e. m ass and viscosity) of the protoplanetary 
disk jKlev et alJ2004h . If a system is a subject to such an adia- 
batic migration, a given value of K result s in unique values of 
the final eccentricities jLee & Pealel2002l) . However, as shown 
above, in the case of HD 12831 1 the eccentricities are varying 
with quite large amplitudes. 

We have performed a set of numerical integrations of the 
general co-planar three-body problem adding non-conservative 
drag forces, varying the value of K. 
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Fig. 1. Top: Time evolution of the resonant angle ©1 obtained 
by numerical integration using initial conditions of Table 1. 
Bottom: The corresponding evolution of the eccentricities. 
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Fig. 2. The behavior of the semi-major axes and eccentricities 
during an adiabatic migration with t„ = 2 x 10 3 years and K = 
5. The migration is stopped between 2 x 10 3 and 3 x 10 3 years 
by applying a linear reduction. 



In these simulations the planets move originally in quasi- 
circular orbits and we start them from a 1 = 4AU and a 2 = 2AU. 
Assuming that only the outer planet is forced to migrate inward, 
we implement a dissipative force which results in a migration 
characterized by the e-folding time r fl = 2 x 10 3 years in the 
semi-major axis, and r e — 4 x 10 2 years in the eccentricity 
of the outer planet (corresponding to K — 5). In order to ob- 
tain the present values of the semi-major axes we slow down 
the migration beginning from t\ = 2 X 10 3 years and decreas- 
ing it linearly to zero until £2 = 3 x 10 3 years. In this way we 
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Fig. 3. The behavior of the eccentricities obtained by a sudden 
stop of the migration of the outer planet. 
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model the smooth dispersal of the protoplanetary disk. After 
this migration process, the system is locked into a deep 2: 1 res- 
onance, the resonant angles ©12, and Am librate around 0° with 
small amplitudes. The eccentricities e\ and e^ are almost con- 
stant (showing only a small amplitude oscillations), where the 
index 1 refers to the inner and the index 2 to the outer planet. 
For the chosen K — 5 we find e\ = 0.46 and <?2 = 0.15, see 
Fig. |2] Smaller/larger values of K yield always systems deep in 
resonance and result in larger/smaller e 12, contradicting Fig.^ 

Clearly, the present behavior of HD 128311 is not the result 
of such an adiabatic migration process alone. In the follow- 
ing, we shall present two additional mechanisms, which may 
be responsible for the large oscillations of the eccentricities and 
breaking the apsidal corotation. 

3.2. Sudden stop of the migration of the outer planet 

Recent Spitzer observations of young stars confirm that the in- 
ner p art of the protop l anetary disk may contain o nly very little 
mass D'Alessi o et al.l (12005): Cal vet et al.l 12005). possibly due 
to photo-evaporation induced by the central star. Thus, upon 
approaching the inner rim of such a disk, the inward migration 
of a planet can be terminated rapidly. 

In order to model this type of scenario, we perform addi- 
tional simulations where the migration of the outer planet has 
been stopped abruptly reaching the actual value of its semi- 
major axis (02 = 1.73AU). We assume that the inner planet 
orbits in the empty region of the disk at a\ = 1.5AU, and we 
start the outer planet from 02 = 4AU forcing to migrate inward 
very fast (j a = 500 years). We find that the present behavior of 
the eccentricities can be obtained by using K — 10. The sudden 
stop of the migration results in a behavior of the eccentrici- 
ties shown in Fig. |3 which is very similar to the observed case 
(Fig.Q. We note that after the above sudden stop of the migra- 
tion the planets remain in apsidal corotation but with increased 
libration amplitudes of the resonant angles. 
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Fig. 4. Time evolution of the eccentricities (Top), resonant an- 
gles ©i (Middle) and Act (Bottom) of the inward migrating 
giant planets, before and after a scattering event with an inner 
low mass planet. 



3.3. Planet-planet scattering event 

The behavior of the eccentricities of the giant planets of the sys- 
tem HD 128311 is ver y similar to that ob served in the system 
around v Andromedae. lFord et al.l J2005h proposed that such a 
behavior in v And is most likely the result of a planet-planet 
scattering ev e nt. In vestigating the resonant system HD 73526, 
iTinnev et alJ J2006I) also suggested that the lack of the ap- 
sidal corotation can be the result of a dynamical scattering 
event. Therefore we investigate whether the present behavior 
of HD 128311 can be modeled by such an effect. We present 
two cases. 

Additional inner planet 

In this case we assume that an additional small mass planet is 
orbiting close to the central star in a quasi circular orbit. When 
the outer giants are far enough from the star the orbit of the 
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small mass planet is stable. However, as the giants migrate in- 
ward approaching their present positions they perturb the orbit 
of the small planet making its motion chaotic, which in long 
term may result in an increase of its eccentricity. Finally, due 
to its high eccentricity, the small planet can suffer a close en- 
counter with one of the giant planets. After the close encounter 
the small planet may be ejected from the system or pushed into 
an orbit with very large semi-major axis. Depending on the ini- 
tial position of the small planet, it can also be captured into a 
mean-motion resonance with the inner giant, and forced to mi- 
grate inward. During the migration its eccentricity will increase 
in a much shorter timescale than by a "pure" chaotic evolution. 
Thus, in this case a close encounter with the inner giant is very 
likely as well. 

In Fig.0]we show the behavior of the eccentricities and the 
resonant angles ©i and Act when the giant planets are migrat- 
ing inwards with initial conditions given in Sect. 13.11 Starting 
the small planet (m = 0.03Mj up ) from a nearly circular orbit 
with a - 0.5 AU, it is captured into a resonance with the inner 
giant. During the migration the small planet's eccentricity in- 
creases, and finally the small planet crosses the orbits of giant 
planets and suffers a close encounter with them. After the en- 
counter ait 2500yrs, the apsidal corotation of the giant plan- 
ets breaks (i.e. ©2 and Act circulate), however the giant planets 
remain in the 2:1 resonance, since ©i still librates around 0°, 
see Fig.|4] 

Additional outer planet 

In this second case we assume that the third planet (m - 
0.03Mj up ) originates from the outer region, and that it ap- 
proaches the giant planets through inward migration where we 
assume that the adiabatic migration of the two giant planets is 
already finished, as displayed in Fig.|2] They orbit for instance 
in a gas-free environment while the small planet is still embed- 
ded in the outer protoplanetary disk. In Fig. [5] we display the 
variation of the eccentricities of the giants after the scattering 
event. The small planet is started from a = 2.6AU and migrates 
inward with an e-folding time r a — 2 x 10 3 years. After scatter- 
ing it is pushed into a very distant orbit a ~ 300AU. We note 
that in this case the giant planets remain in apsidal corotation, 
however with a substantial increase in ©2 and Act. In additional 
simulations we find that apsidal resonance is always preserved, 
however we cannot exclude the possibility of breaking the ap- 
sidal corotation entirely. 



4. Conclusions 

In this paper we investigate evolutionary scenarios, which may 
have led to the observed behavior of the resonant system 
HD 128311, where the two planets are engaged in a 2:1 mean 
motion resonance, but not in apsidal corotation. We assume that 
the two giant planets have been formed far from the central 
star and migrated inwards, due to gravitational interaction with 
the protoplanetary disk. During this differential migration, they 
have been locked into the 2:1 mean motion resonance. At the 
end of their migration we imagine that the system suffers a sud- 
den perturbation. 
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Fig. 5. The behavior of the eccentricities before and after a scat- 
tering with a small mass planet migrating inward from outside 
starting from a = 2.6AU. 

As a first case we study the sudden stop of the migration, 
which may possibly be induced by an inner rim of the disk and 
an empty region inside of it, as indicated by some observations 
of young stars. 

In the second scenario the sudden perturbation is caused 
by a planet-p l anet scattering event, similar to that suggested 
bv iFord et al.l d2005l) in the case of v And. We analyze an en- 
counter with a small (~ 10M ffl ) planet, approaching the two 
massive planets either from outside or inside, which is thrown 
to a large a orbit or directly ejected after the scattering event. 

We find that the sudden perturbation caused by an en- 
counter with an inner small planet is clearly able to break the 
apsidal resonance of the two planets. But also in the other cases 
the orbital behavior of the giants is very similar t o the most re- 
cent orbital data of HD 1283 1 1 JVoet et all2005h . 

The system HD 12831 1 constitutes another example which 
demonstrates the important interplay of migration and scatter- 
ing processes in shaping the dynamics of exoplanetary systems. 
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